Quasiparticle bands and superconductivity in bilayer cuprates. 
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We analyze the generic leatures of the energy spectrum for two coupled Cu02 layers with a 
realistic extended Hubbard model. The quasiparticle bands exhibit flat regions near X(Y) points 
in the Brillouin zone with a large reduction of the bonding-antibonding splitting, and pinning of 
extended van-Hove singularity to the Fermi level, which is more efficient for a bi-layers than for a 
single layer. In contrast to the results with simpler models, the superconducting temperature for 
A^2_y2 pairing is not lowered by the bi-layer hopping. 
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The role of multiple layers in the high-temperature 
superconductors continues to be one of the most in- 
triguing puzzles, with many conflicting proposals and in- 
terpretations. One proposal is that superconductivity 
is enhanced by electronic correlations which are argued 
to greatly reduce single-particle hopping between layers 
while allowing pair-tunneling Jjl| . Angle- resolved photoe- 
mission (ARPES) studies |2|j^ of the energy spectrum 
near the Fermi energy for bi-layer materials can in princi- 
ple answer the questions of the nature of the bonding and 
antibonding bands near the Fermi level. However, one 
recent ARPES experiment |^ resolved only one Cu02 
band at the Fermi surface in BSCO-2212, supporting the 
idea of a greatly reduced interlayer hopping, whereas an- 
other measurement reports two Fermi sheets, in gen- 
eral agreement with the predictions of band structure 
calculations. In addition, there is evidence for important 
effects of coupling between the layers in bilayer systems. 
Neutron scattering experiments see maximum intensity 
for spin scattering for wave vectors the ^ tt/L, where 
L is the interlayer spacing, for both the antiferromagnetic 
insulator and in the superconductor 

A generic feature which has emerged from the ARPES 
experiments |^,|| is a "flat-band" or "extended van-Hove 
singularity" which is "pinned" to the Fermi level for dif- 
ferent hole dopings. Bifurcated saddle points very close 
to the Fermi level had, in fact, been predicted by LDA 
calculations for YBCO-123 Q and caused by dimpling of 
the Cu02 planes (|. The flat band observed by ARPES 
has also been attributed to many-body effects and 
a number of studies have been done for correlated elec- 
tronic states in a single Cu02 plane. The flat region in 
the quasiparticle spectrum has been proposed to be a 
"fcrmion condensate" 0,0, or a non-Fcrmi liquid area 
formed in (q, oj) space near two-dimensional van-Hove 
singularity |l2j. The one-band |^,^ and three-band ||l^] 
Hubbard models as well as t-J model ||^ show a flat quasi- 
particle band just below the Fermi energy, which has an 
"extended" van-Hove singularity near the X(7r, 0)-point, 
due to antiferromagnetic spin fluctuations. The same an- 
tiferromagnetic fluctuations at q « (7r,7r) lead to d^^-y^ 



superconductivity [p^ with a relatively high transition 
temperature. For the two-plane Hubbard model [|l5|jl6[ 
with a simple tight-binding spectrum a reduction of su- 
perconducting correlations was found due to the inter- 
layer coupling. 

In this paper we study the bi-layer Hubbard model 
with realistic LDA derived hopping integrals. Combi- 
nation of this realistic tight-binding (TB) model with 
many-body effects accounts for anomalous properties of 
cuprate superconductors. Thus in the normal state of 
our BSCO model, many body effects strongly reduce the 
splitting between the bonding and antibonding bands in 
the regions in the k-space where the one-body splitting 
is large, leading to flat bands near the chemical potential 
to within an energy of order room temperature. On the 
other hand, in the regions where the many-body effects 
are small, there is also little splitting in the one-particle 
spectrum, due to geometry of the bi-layer bonding in 
the cuprates. We analyze the pinning of the chemical 
potential to these van-Hove singularities in the case of 
a mono- and bi-layer model for different hole dopings. 
A calculation of the superconducting transition temper- 
ature shows robust stability of the A^^-y^ state in the 
"three-dimensional" bi-layer case. 

Let us start with the Hubbard hamiltonian for a bilayer 
Cu02 model: 



H = 



ia.jl3(7 



here tia.jfj are the hopping integrals which gives the en- 
ergy bands en(k), are the site indices inside the plane 
and Oi{[3) are the plane indices, U is the on-site Coulomb 
repulsion and is the chemical potential. Due to the mir- 
ror plane symmetry for the bilayer |^ we can define bond- 
ing (+) and antibonding (— ) electron bands £+(k) and 
£_(k), where k is a vector in the two dimensional Bril- 
louin zone. This symmetry also holds for the quasipar- 
ticles excitation spectrum of the full interacting system. 
For our studies we have used different schemes which 
range from the simplest model (only nearest-neighbor 
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hopping in-plane and between planes) to realistic four- 
and eight-band models, obtained by integrating out the 
high-energy degrees of freedom from the full LDA calcu- 
lations for cuprate superconductors 

For BSCO bi-layer materials we use a four-band Hamil- 
tonian which includes the standard Cu-da;2_y2, and 
0-pa;,Py orbitals, plus the "Cu-s" orbital which has also 
some Cu-d22 character. The last band is needed because 
the standard 3-band model does not give an adequate de- 
scription of the valence band and the Fermi surface [^ . 
The Cu-s orbital provide the 2nd and 3rd nearest neigh- 
bor {t' and t") intra-plane hopping integrals, as well as 
the hoping between the Cu02 planes in the low-energy 
Hamiltonian H. Diagonalization of its first term yields 
the (simplified) "LDA bands" : 

£_,_ (k) — — 2t(cos kx + cos ky) — 4t' cos k^ cos ky — 
2t"{cos2kx + cos2ky) =F t±{coskx — coskyY /A. 

Note that the interlayer hopping is anisotropic with the 
maximum splitting of £+(k) and £-(k) at the X(7r, 0) 
and Y(0, tt) points and no splitting in the direction 
from r(0,0) to M(7r, tt). We chose values of the tight- 
binding parameters for BSCO: t = 0.5 eV, t' /t = —0.3, 
t" /t' = 0.2 and t±_ = 0.15 eV, close to the LDA pa- 
rameters . The total band width is 8i = 4eV, and the 
bonding-antibonding splitting at the X-point is 2ij_ = 0.3 
eV. The screened electron-electron interaction parameter 
for Cud orbitals obtained in the constrained LDA calcu- 
lations for cuprate is 7-8 cV [Q. Since the Cud character 
at the Fermi surface is approximately 65% the effec- 
tive U is reduced by 0.65^ and we used in the present 
study U — 2>.2 eV. This is a bit less than the band width, 
but general results are not very sensitive to the real U- 
values. Note that most many-body Monte-Carlo calcu- 
lations l7|j9|,[l^ use only a simple tight-binding model: 
single - t parameter and isotropic t± . We believe that in 
the cuprate system both many body effects and the band 
dispersion play a crucial role. 

We used the conserving fiuctuation-exchange (FLEX) 
approximation [|9|, which gives an energy spectrum for 
Cu02 plane almost identical to Quantum Monte Carlo 
results . We calculate the self-consistent one-electron 
Green's function in the normal state: 

G'±(k, a;„) ^^^^ _ [e±(k) - /i] - i;±(k, w„) ' 

where I]±(k, a;„) are bonding and antibonding self- 
energies and ojn ~ (2n-|- 1)-kT are the fcrmion Matsubara 
frequencies with n an integer and T the system temper- 
ature. It is convenient to introduce the simple notation: 
k = (k, w„) and = ^ X^k i.^ ' where N is the total 
number of momentum points. The straightforward gen- 
eralization of the FLEX approximation to the case 
of a bilayer two-band Hubbard model gives the following 



equation for the self-energy: 

S±(fc) = Y.^V+{q)G±{k - g) + V-{q)G^{k - q)], 
q 

where the contribution to effective potential V±{q) = 
V^{q) + V^{q) from spin and charge fluctuation are: 

^±-21^77^-^^' ^±-2TT77^-^±' 

here V^'^ = \U'^x± is the double counted contribution 
and V±{q) is defined self-consistently thorough particle- 
hole susceptibility for bi-layer: 

X±{q) - - ^[G+(fc)G±(fc + q) + G-{k)G^{k + q)]. 
k 

We solve the nonlinear integral FLEX-equations using 
the fast Fourier transform (FFT) method |^ on the dis- 
creet mesh of 64x64 momenta in the two dimensional 
Brillouin zone and 600-900 Matsubara frequencies with 
the cutoff of 20-30 eV (which corresponds to the temper- 
ature range of 80-200K). Analytical continuation on the 
real axes was done by Fade approximation. The calcula- 
tions have been carried out for different hole concentra- 
tions: 6 = 0.1 — 0.4 per plane. 
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FIG. 1. The bi-layer spectral function for diflerent direc- 
tions in the Brillouin zone {S =0.25, T=150 K) . 

The spectral function yl(k, lj) = — l/7r/m(G+(k, lj) + 
G_(k, Lo)) for hole doping S = 0.25 per Cu02 and T=150 
K is shown in Fig.|l| for F-X-M directions. One sees that 
there is only a single peak with a "non-Lorentzian" be- 
havior crossing the Fermi level. In fact there are two 
bands with the large renormalization of the interlayer 
splitting from 300 meV to 40 meV, and there is large 
broadening at the actual temperature of the measure- 
ment in the normal state. 
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ceeds that x+('7); for ferromagnetic coupling and near 
q = (tt, tt) the difference could be even an order of mag- 
nitude for t' = t" = 0, reflecting the enhancement of the 
antiferromagnetic fluctuations between the layer . The 
Fermi surface is similar to the LDA, but the large flat- 
ting of the quasiparticle bands and the broadening of the 
spectral function ||T^ is not properly described in LDA. 
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FIG. 2. Energy spectrum for U=0 (TB) and U/0 (QP) for 
the antibonding (top) and bonding (bottom) bands. 
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FIG. 3. Fermi surface for 5 = 0.25, T=150 K: temperature 
broadening (left) and constant occupation number contours 
(right). 



This is illustrated in Fig.^ which shows a compari- 
son of the pure tight-binding spectrum (U=0) and the 
quasiparticle spectrum (QP) for 5 — 0.25 for the bonding 
(b) and antibonding (a) bands together with correspond- 
ing Fermi-surface (zero-energy contour). The quasiparti- 
cle dispersion e2^(k) have been determined through the 
maximum in the spectral function ^±(k, w). Note that 
already the TB (or LDA) bands differ substantially from 
the simple model with only one hopping - t: the position 
of the saddle points at the X and Y points is lower (for 
"t-model" they are located at the middle of the band 
i.e. at 1 eV in Fig.||). This TB-spectrum changes the 
topology of the Fermi surface, and makes the hole doped 
and the electron doped situations quite different (as it 
should be experimentally |^) and suppresses the nest- 
ing near half-filling. The last effect increases the width 
of a susceptibility peak near M(7r, tt) point and reduces 
the antiferromagnetic spin-correlation length from « 3 
to « 1 one lattice spacing. This leads to the suppression 
of the "shadow bands" |21| in the quasiparticle spectrum 
which occur in the simpler model with t' = t" ~ Q and 
smah doping 5 « 0.1 [||. 

A more important consequence of t' and t" is related 
to the increased flatness of the van-Hove singularities in 
the direction towards F, which leads to a larger phase 
space for the spin and charge excitations. This results in 
a drastic renormalization of the energy spectrum (Figj^), 
manly due to spin- fluctuation effects p3| . The quasi- 
particle bands exhibit flat regions near X(Y) points in 
the Brillouin zone with a large reduction of the bonding- 
antibonding splitting. Note that the susceptibility for 
anti-ferromagnetic coupling across the bi-layer x- (?) ex- 



We plot in Fig.||the "temperature" smeared Fermi sur- 
face (/i = ±T for T = 150K) and constant occupation 
number contours: n(k) = ^^^(k, cjn) ||2^ with no—0.5 
corresponding to the Fermi surface. Note the extensive 
regions of nearly equal occupation « 0.5 near the X and 
Y points. It is clear that the difference between the bond- 
ing and antibonding Fermi surface sheets for the BSCO 
model is very small and hard to detect experimentally 
We believe that the "shadow band" obtained in 
the recent ARPES experiment |^ is related to the large 
dark region near X(Y)-point for the antibonding bands. 
This flat region of the quasiparticle spectrum within the 
room temperature scale from the Fermi level gives the 
anomalous linear dependence of the self energy on tem- 
perature and frequency, just from the "phase-space" ar- 
gument |l^j2^ , detected numerically near X and Y points 



In order to show more clearly the formation of the 
flat quasiparticle band in the bi-layer system we plot in 
the Fig.^ the TB and QP spectrum for 5 = 0.35 where 
the Fermi level is approximately between bonding and 
antibonding bands at X. While the Fermi-level crossing- 
points are nearly the same in the TB and QP bands, 
the bonding and antibonding QP bands are "pinned" to 
the Fermi level within room temperature scale. In other 
words, the system has increased the spin and charge fluc- 
tuations by pinning the saddle points to the Fermi level 
and forming the extended van- Hove singularities instead 
of the standard renormalization of the spectrum only in 
the small region near the Fermi surface. 
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FIG. 4. Energy spectrum for TB and QP near X-point. 
Arrows indicate the room-temperature scale, r/2 = (7r/2,0), 
M/2 = (7r,7r/2) (5 =0.35, T=150 K). 

We have analyzed the effect of the saddle-points pin- 
ning to chemical potential for the monolayer {t±^ = 0) and 
the bilayer models. The energy position of the saddle- 
point at X with respect to the Fermi level for the TB- 
and QP-spectrum as functions of the hole doping shown 
in the Fig.|[ The antibonding QP-band for the bi-layer 
model (Bi-a) stays within the ±5 meV for large hole con- 
centration range 0.1 < 6 < 0.3 per Cu02- Such a pinning 
effect is less pronounced in the monolayer (Fig.^). 




FIG. 5. Position of the saddle point at X relative to Fermi 
level for the TB and QP bands. Note the different scales. 

The many-body reduction of the interlayer splitting 
also increases the quasiparticle density of states in the bi- 
layer and could help the formation of "three-dimensional" 
superconducting state. We find the superconducting 
transition temperature by solving the linearized Eliash- 
berg equation for the singlet order parameters in the bi- 
layer model I p^ . The critical temperature is defined via 
the maximum eigenvalue: Amax(7c) > 1 of the following 
superconducting kernel: 



Rk'. 



Vl''G+{k)G+{-k) Vl''G+{k)G+{-k) 
Vl''G-{k)G-{-k) Vl''G-{k)G-{-k) 



where ^1^= = F^fc' - fc) - (fc' - fc) -f L/. We find that the 
superconducting d^2_y2 state, which is symmetrical with 
respect to the mirror plane between the bi-layer, is indeed 
quite stable. The transition temperature for 5 = 0.2 is 
about 90 K (or Tc=0.015t), which is almost the same as 



for single layer (within one percent). The theoretical su- 
perconducting temperature for a LSCO-model |2|] with 
a nested Fermi surface is a bit higher. If one takes into 
account the large order-parameter fluctuations which are 
expected to reduce Tc in the single layer we can 
expect higher Tc's in the bi-layer systems. 

We have repeated the calculations for the YBCO eight- 
band TB-model |^ and found that the "bifurcated" sad- 
dle points in the TB-spectrum helps to create the ex- 
tended anisotropic van-Hove singularity. In the YBCO 
model interlayer splitting is reduced by only a factor of 
two and could be in principle detected by ARPES mea- 
surements. 

In conclusion, we have shown that the generic feature 
of bi-layer as well as mono-layer cuprate is the forma- 
tion of extended van-Hove singularities near the chemical 
potential on the scale of the temperature, which occurs 
more readily in bi-layers. This leads to large areas of the 
Brilluoin zone with nearly constant occupation « 0.5, 
which can explain the anomalous normal state proper- 
ties and lead to enhanced dx2-y2 superconductivity. 
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